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Abstract 

Mueller's dipole formulation of onium-onium scattering is used to study unitarity correc- 
tions to the BFKL power growth at high energies. After a short discussion of the spatial 
distribution of colour dipoles in a heavy quarkonium and the associated fluctuations, re- 
sults are presented showing that the one and two-pomeron contributions to the total cross 
section are the same at a rapidity Y ~ 14. Above this rapidity the large fluctuations in 
the onium wave function cause the multiple pomeron series to diverge. Resumming the 
series allows one to show that unitarity corrections set in gradually for the total cross 
section, which is dominated by rare, large, configurations of the onia. The elastic cross 
section comes mostly from much smaller impact parameters and has significant unitarity 
corrections starting at a rapidity y ~ 8. 
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1 Introduction 



The BFKL pomeron [||, |2|, causes the ampUtudes for certain high energy processes, which 
are also hard processes, to rise as a power of the centre of mass energy. At very high energies 
this rapid rise will cause the amplitudes to violate unitarity. Unitarity is the statement that 
for any given initial states and impact parameter, the probability of interaction cannot be 
larger than one. Because it is a limit on the interaction for each impact parameter, one cannot 
tell when an amplitude violates unitarity without knowing how the amplitude is distributed 
in impact parameter. This has also been one of the main uncertainties in calculations (for 
example [^) of the shadowing corrections to BFKL evolution using the GLR equation Q. 
However knowledge of the average interaction for a given impact parameter is insufficient. The 
wave functions for each of the colliding objects consist of a variety of states each with some 
amplitude. What can happen is that for a given impact parameter the average interaction 
can be dominated by occasional rare combinations of states which have a large interaction 
and therefore large unitarity corrections. However because these configurations are rare, the 
average interaction can still be small and one can be misled into thinking that there should 
be no unitarity corrections. 

One of the simplest frameworks in which to study the BFKL pomeron, and the onset of 
unitarity corrections, is the dipole formulation, proposed by Mueller 0, ^, ^. For sufficiently 
heavy onia, as, the strong coupling constant, is small and one can apply perturbation theory. 
One starts with a qq state and considers the correction to the wave function which comes 
from adding one gluon, causing the initial qq colour dipole to branch into two colour dipoles 
{qg and gq). In the small-x, leading Nq approximation, each of these dipoles can then branch 
again, independently, and the process repeats itself, building up a cascade of colour dipoles. 
One can then calculate the onium-onium elastic amplitude from the interactions between the 
dipoles of the two onia. The lowest order interaction is the exchange of a colour neutral pair 
of gluons. This is one pomeron exchange. Exchange of larger numbers of pairs of gluons 
corresponds to multiple pomeron exchange. This approach is particularly suited to the study 
of unitarity because it is formulated in impact parameter and also because it contains all the 
information on the fluctuations in the onium configurations. Further, because the onium wave 
functions are each evolved to only half the total rapidity, unitarity corrections to the cross 
section should set in long before the density of partons rises to the point where corrections 
to the wave functions need also be taken into account, i.e. long before the wave functions 
saturate. 

In this paper, a Monte Carlo simulation will be used to calculate the full, unitarised, 
amplitude for high energy onium-onium elastic scattering. Section § gives a brief review 
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of the dipole formalism. Section ^ discusses certain properties of the onium wave function, 
showing that the distribution of dipoles is more central than had been suggested in |^], and 
giving more information on the exponential tails of the fluctuations in dipole density, which 
were first mentioned in [10|. These fluctuations have as a consequence that the /c-pomeron 
series is divergent, though it can be resummed by interchanging the order of the summations 
over k and over configurations of the onia [Q]. This and other consequences of fluctuations 
are illustrated in section ^, where the main results of the paper are presented: the fc-pomeron 
and unitarised amplitudes for onium-onium scattering, as a function of rapidity. One and two 
pomeron contributions to the total cross section are seen to be equal at a rapidity of 14, but 
the unitarised total cross section carries on growing well beyond this point. The elastic cross 
section has much stronger unitarity corrections and the behaviour of the power of the growth 
of the unitarised and 1-pomeron elastic cross sections differ significantly above 1" = 8. 



2 Onium-onium scattering 
2.1 Single pomeron exchange 

The amplitude for exchange of a single pomeron between two onia moving fast in opposite 
directions, with a centre-of-mass energy is: 

^(r, Y) = -iJ d^bd^b' J dzdz'$(b, z)$(b', z')F{h, b', r, Y) (1) 

The rapidity Y is approximately log s/M^, with M being the onium mass. The amplitude 
is determined for a fixed relative impact parameter r between the onia and ^(h,z) is the 
square of the heavy quark-antiquark part of the onium wave functions with b the transverse 
separation between the quark and antiquark and z the longitudinal momentum fraction of 
the antiquark. For single pomeron exchange, F = F^^\ with 

(b, b', r,Y) = -J ^^d^Kd^n'fiK - R', c, c')n(c, 6, R, y)n{d , 6MR' - r| , y - y) 

(2) 

The factors n(c, 6, i?, y) are the density of dipoles of size c at a distance R from a parent onium 
of size h after evolution through a rapidity y. The terms 1 /2'K(? are just to compensate for the 
normalisation of each wave function (which from now on refers to the gluon distribution of 
an onium of fixed size). Note that for one pomeron exchange the division of rapidity between 
the two onia (the choice of y) is arbitrary. The amplitude (divided by (— «)) for interaction 
between a dipole of size c and one of size c' whose centres are separated by r, /(r, c,c'), is 
the following: 
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dV-/(r,c,c') = f ^ 



AH 



g-i£.c'/2 _ gi^.c'/2 



gi(q-£).c/2 _ g-i(q-£).c/2 



^_i(q_^).c'/2 _ gi(q-^).c'/2 



(3) 



This is considered in more detail in the next section. The normalisation is chosen such that 
the total cross section is: 



atot{Y) = 2 j d?TluiA{T,Y) 



(4) 



while for the elastic cross section 



d(7eZ _ J_ 
di ~ 47r 



(5) 



with t = —q^ and 



A{q,Y) = j d\e^'^-''A{v,Y). (6) 

The term total amplitude will be used to refer to F{b,b' ,r,Y) integrated over r, which is 
equal to half the total cross section (ignoring the integral of eq. (|^) ) . 

In the limit r ^ b, b' , and with logr^/66' <C kY, one has the following expression for F(i) 



r'^ \ bb' I [-KkyYI^ 

where (a-p — 1) = 41og 2a5'A''c/7r and k = lAasNcQi^j/T^, C{f^) being the Riemann zeta 
function. 

It is convenient to work in impact parameter and with this normalisation, because the 
unitarity condition is then simply |-F(r, y)| < 1. It is the exponential dependence on Y which 
causes this condition to be violated as Y increases beyond a certain limit. 

2.2 Dipole-Dipole interaction 

In the dipole-dipole interaction, /(r, c, c') was approximated by: 

/(r, c, c') = <52(r)7r44 {l + log ^) (8) 

This is correct if one averages over angles and integrates over r, and is therefore perfectly 
adequate for studying single pomeron exchange. However to study multiple pomeron exchange 
one needs to know the full expression. It is obtained by performing the inverse Fourier 
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transform in eq. (^), noting that f{q), the Fourier transform of the potential, is a convolution 
of two identical terms, each of which gives the single gluon potential between a pair of dipoles. 
The result is that /(r) is then simply the square of the two-dimensional dipole-dipole potential: 



/(r,c,c') 



a 



2 r 

s 



log 



|r + c/2 - c'/2\ |r - c/2 + c' /2\ 



(9) 



|r + c/2 + cV2||r - c/2 - cV2| 
as would be expected for the exchange of two gluons between a pair of dipoles in two dimen- 
sions. 

2.3 Multiple pomeron exchange 

To recover the unitarity condition it is necessary to include higher order terms. The ones 
which will be considered here are /c-pomeron exchange. This corresponds to the exchange of 
k colour neutral pairs of gluons between the two onia. As is discussed in [§] and reviewed here 



in section 2.4, at rapidities where unitarity effects are setting in, it should be safe to neglect 
corrections to the wave functions which would arise due to saturation of the wave functions. 
To evaluate multiple pomeron exchange, it is no longer adequate to use the average dipole 
distribution: it is necessary to know the higher moments of the wave functions 

In , the equations governing the evolution of the wave function are expressed using an 
operator formalism: the operator a^(r,b) ((i^(r,b)) creates a dipole at position r and with 
size b, in the right (left) moving onium. The basic vertex of the dipole picture is then Vi: 



Vi = ^J ^d2boid2b2dW(r + ^,bo2)at(r + ^,bi2)a(r,boi). (10) 



^2 = -%2^ / ^2^d2boid2b2d2ra(r,boi)at(r,boi) (11) 



"02"12 

This expresses the branching of a colour dipole of size bgi at position r into the two dipoles 
which are formed by the production of a gluon at b2. Together with this there are virtual 
corrections 

asNc f _6g 
I, 

^2*^12 

which ensure the conservation of probability: when a dipole branches into two, the original 
dipole is lost. For a given dipole, the branching process and the virtual corrections depend 
only on its size, not on its position or on the other dipoles that are present. The probability 
of finding a configuration made up of dipoles at positions and sizes {r 1 , Ci ; . . . ; r^, c^} after 
the evolution through a rapidity y of an initial dipole at tq of size b = boi is then: 

P({ri, ci; . . . ; r„, c„}, Fq, b, y) = (0|a(ri, Ci) . . . a(r„, c„)e^^«at(ro, b)|0) (12) 
where Vr = Vi + V2. There is a similar expression for the left-moving onium. 
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The amplitude F^^\r,Y) for the exchange of k pomerons is then, analogously to eq. |2|, 

FW(r,y) = (0|e'^i+'^i-^^VT^ef^^+(^-2')^«dt(r + ro,bL)at(ro,bK)|0). (13) 
The term ai is 

ai= J d^rd^c a(r,c) (14) 

with an analogous expression for di. The operator / gives the sum of all the dipole-dipole 
interactions occurring for that impact parameter: 

f = J d^Rd^R'd^cd^c'f(R-R',c,c')d\K,c)d(R,c)a'<(R',c')a{R',c'). (15) 

Summing together all orders of pomeron exchange, one obtains an expression which explicitly 
satisfies the unitarity bound: the 5-matrix is 



S{h,h',r,Y) = 1 + F(b,b',r,y) = (0|e'^i+'^ie-^e^^^+(^-^)^«d^(r + ro, b)t(ro, b')|0). (16) 

Where it would otherwise be unclear which quantity is being discussed, the unitarised am- 
plitude F will be referred to as This can be cast into a form more suitable for use 
with a Monte Carlo simulation. Let 7 be a particular dipole configuration for the left moving 
onium, which contains dipoles of position and size (ri, Ci) . . . (r„^ , c„^). For an onium at 
impact parameter r, of size b, and evolved to a rapidity of Y, the probability of finding such 
a configuration is defined as P^(h,r,Y). In these terms, the 5'-matrix becomes: 

S{Y, r, b, b') = ^ P^(ro, b, y)Py (ro + r, b', Y - y) exp(-/^,y ). (17) 

7.7' 

The definition of f^^y is: 

/7,7' = EE/(r-r;,c„c;.) (is) 
i=i j=i 

With this form for the interaction, the cross section can be obtained by randomly producing 
the configurations for each onium, and working out f^^y for each point in impact parameter. 
This is repeated with many configurations to obtain the average S-matrix. Since f^^y is 
always positive, it is obvious from eq. ^ that the unitarity bound is satisfied. 



5 



2.4 Division of the rapidity interval 

As already mentioned, for one pomeron exchange, the division of the rapidity interval between 
the two onia is arbitrary. The onium-onium scattering amplitude is independent of y in eq. (|2|). 
This is because y is not a physical parameter. For multi-pomeron exchange, though y is still 
not a physical parameter, it does become a parameter of the approximation, because one is 
ignoring corrections to wave functions which have been evolved to y and {Y — y) the 
fractional correction to the forward amplitude from 2 pomeron exchange is ©(a^e^"'''"^-*^). 
For any given part of the wave function, the leading correction to that part of the wave 
function is of the order of 0(a|e^'^'''~"'^''^) (or 0(a|e("'''~^)*^^~^)) for the other wave function). 
If one divides the rapidity equally between the two onia, y = corrections to the wave 
functions do not become important, in principle, until twice the rapidity where corrections 
appear in scattering. However if one divides the rapidity interval asymmetrically (e.g. y ~ y), 
corrections to the wave function which have been neglected are of the same same order as 
corrections to scattering amplitude and one will underestimate the total unitarity corrections. 
This of course ignores any of the details of the impact parameter distribution of the wave 
function, as well as fluctuations which might cause occasional larger dipole densities, but the 
basic argument should still be valid. 



3 Some properties of the wave function 
3.1 Monte Carlo calculation of the wave functions 

The first stage in a Monte Carlo calculation of onium-onium scattering is to simulate the 
branching processes which produce the wave functions of the two onia. For each onium, 
one starts with an initial dipole of size b and produces a gluon with a random position 
whose distribution is determined by eq. ([To|). The rapidity of the gluon is chosen with an 
exponentially decaying distribution in accordance with the virtual corrections of eq. (11). The 



procedure is then repeated with the two new dipoles, and then the dipoles they produce, etc., 
until the rapidity of any new gluons would exceed the rapidity to which one is evolving. The 
main difficulty in implementing this procedure is that the integrals in eqs. ( p!o|) and (^) are 
divergent for small dipole sizes. This divergence has to be regulated, and the most appropriate 
method for a Monte Carlo simulation is the introduction of a lower cutoff (qow) on the dipole 
size. In any observable quantity the divergence cancels out between the real and virtual terms, 
and the dependence on the cutoff should disappear as it is taken to zero. For a finite cutoff, 
though, any observable will be somewhat modified. The problem with reducing the cutoff 
to zero is that the number of dipoles in the wave function is oc \/ciow- When considering 
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onium-onium collisions the time taken to calculate an amplitude is the proportional to the 
product of the number of dipoles in each onium, and so is oc As a compromise, all 

evolutions will be carried out with a lower cutoff of O.Olfe. This will give good results over 
most of the rapidity range. The value of as used will be 0.18, corresponding to a scale of 
about lOGeV. 

Before looking at onium-onium scattering, there are certain properties of the wave function 
which are worth highlighting. 

3.2 The spatial distribution of dipoles 

An approximate form for the average dipole distribution has been given in [^: 




Figure 1: The spatial distribution of dipoles in the onium wave function, for y = 14 with 
c = b. 

This is the number density of dipoles of size c, at a distance r from the parent dipole of size 
b, after evolution through a rapidity y. The equation is valid for r ^ b,c and log r"^ /be <C ky. 
One of its main characteristics is that the maximum of the dipole distribution is at large r 
(see the dashed line of figure |l|) . This would imply that scattering is taking place at values of 
r where the average dipole distribution would be quite dilute (though fluctuations would still 
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have to be taken into account). It would also have the consequence that one might reach the 
dangerous infra-red region, where non-perturbative effects come into play, somewhat sooner 
than expected. The Monte Carlo results (the + symbols on figure ||), however, show that the 
distribution is far more central. 

The source of this discrepancy is a missing portion of one of the integrals used in the 
derivation of eq. (p!9[). The details are discussed in Appendix A, where it is shown that a 
better approximation has an extra factor in the logarithms: 

nic h r v) -^loJ^-^] exp[(ap-l)^-log^(16rV6c)/%] 

This modification also applies to the one-pomeron exchange amplitude, which is why eq. (|^ 
differs from the result derived in 

3.3 Localised fluctuations 

As already discussed, it is not enough to know just the average distribution of dipoles within 
the wave function. When calculating multiple pomeron exchange it is necessary to understand 
the fluctuations. In |1C], it was found that there are very signiflcant {P{n) oc exp(— log^n)) 



fluctuations in the total number of dipoles of a given size in the onium wave function. The 
source of these fluctuations is the occasional production of a large dipole which cascades 
down to give many smaller ones. For multi-pomeron corrections to onium-onium scattering, 
however, these are not very important, because the dipoles produced from the large dipole will 
be spread over a large region, giving quite a dilute dipole density. But it was noted that the 
number of dipoles in a localised region also has significant fiuctuations, this time exponential. 
As is discussed in ||8| (and as will be reviewed in section |3.4D an exponential distribution in the 
dipole density can cause the series for the amplitudes of /c-pomeron exchange to be divergent. 
Therefore these exponential fluctuations warrant more study. 

A detailed derivation of the exponential nature of the fluctuations is presented in appendix 
B. The basis of the derivation is that at very large Y, the c and b dependences factorise, as 
can be seen from eq. (pO|). The q''^ moment of the distribution of dipoles of size c inside a 
region of size p, from a parent of size b which is also in the same region, can then be expressed 
in the form: 

n(^\c,b,y)=egqlA'i{c/p) (^^^^^("^-Ds/, (21) 

where, for the purposes of the analysis, A can be any function, though it will actually be 
proportional to p/c. All the moments are defined by the coefficients eg and the exponents 
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Uq which are approximately independent of h. Note that this form exphcitly exhibits KNO 
scahng [11| in c and y. One then obtains the following relations, valid for large q: 



2glog2 

and 



^9 - 2 - TTTTTTTT' (22) 



11^ 1 



e. 



^-41og2g^'^''^-^, + z.,_,-2- 
Therefore for large g, one has Cg cx C"^, where the exact value of C depends strongly on the 
behaviour of Cq for low q. This gives for the probability of finding n dipoles of size c within 
a region of size p (in the limit of large n): 



Pn{c,b,p,y) oc exp 



n 



(24) 



A{c/p)Ce(''^-^^y_ 

This agrees with the Monte Carlo results to within the uncertainties, which are due to the 
accessible values of y being non-asymptotic and also because the above formulas are valid 
only for large q, while the value of C retains a dependence on the behaviour of Cq at smaller 

q- 

Note that the shape of the distribution for large n is independent of the size b of the initial 
dipole: i.e. the shape of the tail of the distribution is independent of the starting conditions 
of the evolution. For example it is possible to produce high density fluctuations anywhere 
in impact parameter, as long as a dipole of size ~ p is produced at that point in impact 
parameter, early on in the evolution, so that it has a long range in y to produce many child 
dipoles. The only aspect of the tail of the distribution which does depend on the starting 
conditions is the normalisation: in the above example, the probability of producing a dipole 
of the appropriate size, at the correct impact parameter, depends on the position and size of 
the parent dipole. 

3.4 Consequences of fluctuations: a toy model calculation 

The significance of the exponential density fluctuations is that in Mueller found that 
exponential fluctuations in a toy model without transverse dimensions led to divergence of 
the series for the amplitude of A:-pomeron exchange. Consider an onium at a rapidity Y/2, 
which has a mean number of dipoles p. (= e^"''"^)^/^), with an exponential distribution for 
the probability of there being n dipoles: 

P„(y/2) ^ ie-"/'^ (25) 
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The amplitude F^^'^ for the interaction of two onia by exchange of k pomerons is 



{—a^mn)^ 



(26) 



kl 



n 



This is vahd for k <^ 171,11. Substituting the exponential probability distributions, and ap- 
proximating the sum by an integral gives 



which for a^/i > 1 is divergent. The same idea applies in 4 dimensions, as one may see 
by restricting consideration to the interaction of dipoles of size c ~ 6 within regions of size 
b in each onium (this forces all the interactions to take place at roughly the same impact 
parameter, so that they all contribute to multiple pomeron exchange). These restrictions are 
appropriate at small impact parameters where the interaction will be dominated by the same 
region in each onium. 

4 Monte Carlo study of onium-onium scattering. 

Once one has used a Monte Carlo program to generate a pair onium configurations, as de- 
scribed in section ^, one can go on to look at the interaction between them. This is illustrated 
in figure ||. The top two plots show typical dipole configurations of two onia. The positions 
of the initial quark and antiquark are represented by black discs and each 'cigar' shape rep- 
resents a colour dipole. The evolution in rapidity results in a long chain of dipoles stretching 
between the quark and anti-quark. The lower left hand plot shows the one-pomeron exchange 
amplitude as a function of impact parameter (the position of the second onium relative to 
the first). One can see that the largest interactions occur at an impact parameters where 
large numbers of dipoles overlap. At impact parameters where no dipoles overlap, there is no 
interaction. The plot is for a high rapidity (Y = 20) and most of the interaction is above 1, 
and so violates the unitarity limit. The lower right hand plot shows the unitarised amplitude, 
which is quite flat and close to 1: as one would expect, at large rapidities, the interaction is 
mostly black. 

When determining average amplitudes, two methods have been used. One, as in figure 
involves determining the interaction at each point of a grid, whose size is chosen to encom- 
pass all of the interaction. This gives an estimate of the total interaction for each Monte 
Carlo event. The second approach is to determine the interaction only along a radial line 
in impact parameter, from which, after averaging over many events, one can reconstruct the 
total amplitude. Because one is sampling the onium-onium interaction at fewer points in this 




(27) 
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one pomeron amplitude unitarised amplitude 

Figure 2: The dipole structure of two onia, each evolved up to rapidity y = 10 (each 'cigar' 
shape represents a dipole). The black discs represent the initial quark and anti-quark. The 
bottom plots show the amplitude for their interaction as a function of relative impact param- 
eter. 
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second method, the computational time per pair of configurations is substantially reduced 
and one can look at a much larger number of configurations. Full details of the Monte Carlo 
procedure will be presented elsewhere [ 12 1. 



4.1 1 pomeron, 2 pomeron and unitarised amplitudes 

Before drawing any conclusions from the Monte Carlo calculations, it is necessary to check 
that the answers obtained are in good accord with expectations. All the results which will be 
given here will be for the collision of onia of the same size, b. 




Y 

Figure 3: The total amplitude for onium-onium scattering as a function of rapidity. See text 
for details. 

Figure ^ shows Monte Carlo results for 1 pomeron exchange, compared with a result from 
the numerical solution (as in |l^) of the differential equation governing the evolution Q. 
Slight mismatches between Monte Carlo line segments are visible at y = 12 and Y = 16: the 
points between Y = 12 and 16 all share the same evolution up to y = 12, so that the errors 
due to fluctuations in the onium wave function are correlated, making it easier to see the 
trends in the evolution. The points from y = 16 to 20 share their evolutions up to y = 16, 
but they use a different set of evolutions from the points y = 12 — 16, hence the mismatch 
at y = 16 (and analogously at y = 12), which is a measure of the statistical error. Note 
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that the fluctuations in the wave functions increase at larger Y, while the accessible statistics 
go down (because of the larger number of dipoles in the onia), so that the statistical error 
increases at large Y. This is translated into a downwards shift of the amplitude, because one 
usually misses out on fluctuations with very large amplitudes which would tend to increase 
the average amplitude. In addition, at higher rapidities a suppression of the same order sets in 
due to the lower cutoff on dipole size. These small factors aside, the Monte Carlo simulation 
gives a good estimate of the amplitude. 

The main point of flgure ^ is to show the relation between the two-pomeron and the one- 
pomeron amplitudes. This is of interest because the rapidity at which the two amplitudes are 
the same should be related to the rapidity at which unitarity corrections become important. 
The cross over occurs at y ~ 14, corresponding to an energy which is well beyond the reach 
of current accelerators. 

In [^, Mueller and Patel derived the following relation for the total two pomeron ampli- 
tude: 

The expression which was given for A was left unevaluated. Eq. (28) is plotted in flgure ^, using 
the Monte Carlo value for F^^\ with a fltted value of A = 6.9 x 10^ (including uncertainties 
due to the cutoff, there is an error on this of about 5%). The agreement is very good, except 
at low rapidities where non-asymptotic contributions are important. 



The results of sections and |3^ suggest that, beyond the point where one and two 
pomeron exchange cross over, the /c-pomeron exchange series will be divergent and hence of 
no use in calculating the total cross section. This is clearly visible in flgure Q which shows 
1 to 4-pomeron exchange for Y = 12-16. Above Y = 14 the series looks divergent, while at 
lower rapidities, it looks asymptotic, with 2 and 3 pomeron exchange crossing over before 1 
and 2 pomeron exchange, and so on and so forth for higher k. 

One can test the ideas of section 3.4 in more detail by plotting F^'^^r = 0)/F^''-'^\r = 0). 



In the toy model discussed there, this would be equal to k^asfJ-)^, giving a straight line 
when plotted against k was the mean number of dipoles in the onium without transverse 
dimensions). The same should apply for the case with 2 transverse dimensions, except that 
{as^j)"^ will be replaced by some factor which should be proportional to F*^^^(r = 0). The 
complications due to the transverse dimensions mean that it is not possible to make any 
prediction about the constant of proportionality. 

Figure |5| shows the ratio F^'^\r = 0) / F^'^^^^r = 0) for Y = W. A low value of Y is 
chosen to allow the greater statistics which are needed to sample the tails of the probability 
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Figure 4: Forward, fc-pomeron amplitude for onium-onium scattering as a function of rapidity, 
as obtained from the Monte Carlo simulation. 



distributions, which are relevant for large k. As one can see, the ratio has a clear linear 
dependence on k. The intercept with the k axis is at non-zero k, indicating that the k 
pomeron amplitude behaves as {k + n)!, rather than k\, where —n ~ —3 is the value of the 
intercept. This difference means that, after taking into account the effects of the transverse 
dimensions, the fluctuations of the wave functions effectively deviate from exponential at 
low multiplicities, which is not unreasonable (the calculation showing the fluctuations to be 
exponential was valid only for the tails of the fluctuations). If one examines plots of the same 
ratio for other values of Y (in the range y = 8 to 12) one finds that the slope of the straight line 
is approximately equal to (3F^^^ (r), where /? ~ 0.2. From the linear dependence of the ratio on 
k, and the behaviour of its slope with Y, one concludes therefore that the fc-pomeron exchange 
amplitude at r = is roughly proportional to {k + n)\[PF^^\0)]'' , confirming that when 
unitarisation corrections start to become important, the multiple pomeron series diverges. 

So at high rapidities, to understand unitarity corrections it is necessary to use the full 
unitarised amplitude, resumming multiple pomeron exchange before summing over configu- 
rations of the two onia. The results of this are shown in figure ^. Unitarity corrections set 
in very gradually: even when the amplitude for two pomeron exchange is many times larger 



14 



^ 1-1 




r ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■! 

1 2 3 4 5 6 7 8 

k 



Figure 5: Ratio of k to {k — l)-pomeron exchange amplitude for impact parameter r = 0. 



than the one pomeron amplitude, the unitarised amplitude carries on growing relatively fast. 



At a rapidity of y ~ 19, the reduction in the effective power is only about 25% (see figure 10). 
The key to understanding why unitarity corrections set in so gradually is in the the profile of 
the amplitude in impact parameter. 

4.2 The profile of the amplitude: F(r). 

The profile in impact parameter of the amplitude, F{r), allows one to determine which region 
contributes most to the cross section, and also which region which has the largest unitarity 
corrections. The top plot of figure shows F{r) against r. As one would expect, the amplitude 
is largest at small r and dies off quickly at larger r. However, to understand the contribution 
to the total amplitude, it is necessary to weight the amplitude with integration area: 

F = J d\ogr 2TTr^F{r). (29) 

The lower plot of figure |^ shows the amplitude weighted with , so that the total amplitude 
is proportional to the area under the curve. The largest contribution to the total amplitude 
does not come from the region where F{r) is largest: the enhancement means that radii 
of the order of 2 or 3 times the parent size are the most important for the total amplitude, 
even though the average amplitude is considerably smaller there than at r = 0. 

The two-pomeron amplitude is also shown in figure 0. With increasing r, it dies off faster 
than the one pomeron amplitude, which means that when weighted with r^, it is more central. 
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Y 



Figure 6: Total 1-pomeron and unitarised onium-onium scattering amplitudes, as obtained 
from the Monte Carlo simulation. 



This is as suggested by Mueller in and the consequence is that the one-pomeron and 
two-pomeron amplitudes are dominated by different regions of impact parameter: the two- 
pomeron amplitude can become large without there being large unitarity corrections, because 
the two-pomeron contribution (and unitarity corrections) grows mostly near r = 0, while the 
total amplitude comes predominantly from larger radii. Note that while F^'^\r) / F^-^\r) 
decreases with r, the ratio F^'^\r) / [F^^^r)]"^ increases quite rapidly with r, as was predicted 
in Q (in fact F^'^\r) decreases even more slowly with increasing r than was suggested there). 

Looking again at figure ^, one sees that the largest unitarisation corrections are at also 
small r, so that the unitarised amplitude, once weighted with r^, has its maximum at larger 
r than the one pomeron amplitude. However, the unitarised amplitude is still considerably 
suppressed (compared to the one pomeron amplitude), even at radii where F^^\r) <C 1. This 
is indicative of large fluctuations, or alternatively that eikonalisation is a poor approximation 
(and is also related to the strong rise in the ratio F^^\r)/[F^'^\r)]'^ with increasing r). The 
eikonal approximation is valid in the limit of no fluctuations: if F^^^ (r) is the same for each 
event then unitarisation can be done after averaging over events: 

_ punit^^-^ ^ -F'^'^r) = 1 - e-^'''(^\ (30) 

bearing in mind that F^^^ is negative. Figure I shows the ratios F"™7f(i) and 

For very small r where the amplitude is quite black, the eikonalised and unitarised amplitudes 
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Figure 7: The top plot shows the average amplitude as a function of r. In the lower plot, the 
amplitude is weighted with so that the area under the curve is proportional to the cross 
section. Y = lA. 



are both strongly reduced, and to a similar extent. They do differ slightly, and this is the form 
of failure of eikonalisation found in the toy model . However if one examines the radii which 
contribute most to the total amplitude (r ~ 26 — > 36), eikonalisation would predict only very 
small corrections, whereas the unitarised cross section is reduced by a factor of two roughly, 
compared to the one pomeron value. This means that the configurations contributing most 
to the one-pomeron cross section actually have relatively large amplitudes. Therefore for the 
average amplitude to still be low, the one-pomeron amplitude must be dominated by rare 
configurations. 

The lower plot of figure ^ shows that this is in fact the case. The solid curve is the 
probability distribution of the amplitude F^^\r) for |r| = 26. The dashed curve shows 
this probability distribution weighted by F^^\r), normalised so that the total area under 
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Figure 8: Top plot: the ratio of the unitarised and eikonalised to the one pomeron ampUtude, 
as a function of r. Bottom plot: the probability distribution of F^^^ for a fixed impact 
parameter, |r| = 26. Both plots are for Y = 14 



it is 1. The area under a particular region of the curve is indicative of how much that 
region contributes to the average amplitude at r. The unweighted distribution is very sharply 
peaked at F^^) = 0. Yet the weighted distribution is spread over a wide range of F^^^ : 80% 
of the amplitude comes from from only 10% of the configurations. This is at the value of r 
which contributes most to the total amplitude. There is also a contrast between the curves 
for 1 pomeron and 2 pomeron (i.e. weighted with [F^^^]'^(r)) exchange, however it is less 
marked: the region contributing 80% to the two pomeron amplitude contributes 40% to the 
one pomeron amplitude, though this region comes from only 2.5% of configurations, so that 
the two-pomeron amplitude here is dominated by somewhat rarer configurations than the 
one-pomeron amplitude. The same phenomenon, of rare configurations dominating one and 
multi-pomeron exchange is even more accentuated at larger r. It is also worth pointing out 
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that since the total amphtude comes mostly from moderate values of r, it too is dominated 
by large fluctuations: at y = 14, two thirds of the total one-pomeron amplitude comes from 
only 10% of events. 

The main reason why the amplitude at a given impact parameter is dominated by such 
a small fraction of configurations is that for an impact parameter somewhat away from the 
onium, dipoles are most likely to be first produced there later on in the evolution. But if they 
appear there only late then they won't have a sufficient rapidity range to build up in numbers. 
However if a dipole is produced there early on (an unlikely occurrence) then it will cascade into 
many dipoles and give a large interaction. These rare, (and at larger Y, strongly unitarised) 
configurations dominate the amplitude, so that even though the amplitude will on average 
be small, it can have significant unitarity corrections. In addition these rare configurations 
will not have circular symmetry, so that even when a state of size r is produced, a point r 
will only have an interaction if the orientations of the configurations are appropriate, further 
increasing the fluctuations in the amplitude at a given point. 

What allows the amplitude to continue growing at large rapidities, is that, for any given 
impact parameter, the probability of the configurations with large interactions can carry 
on increasing until they become common. But there will always be larger radii at which 
large interactions are not yet common, so that the total amplitude can carry on increasing 
through an effective increase in the total area of interaction. This ignores of course the effect 
of saturation of the wave function, or confinement effects placing a maximum limit on the 
dipole size. 

4.3 Elastic scattering 

The unitarity corrections to the total amplitude (and correspondingly, total cross section) 
set in quite slowly, because the increase in the effective area of interaction allows the total 
cross section to carry on rising. However the amplitude is strongly unitarised at small impact 
parameters, and concentrating on this region should make it much easier to see the onset of 
unitarity. A crude way of doing this is to look at the elastic onium-onium cross section. An 
expression for the differential elastic cross section was given in eq. (^), but it is best to look 
initially at the integrated elastic cross section which can be written: 



One should remember that A{r, Y) is obtained from -F(b, b', r, Y) by integrating over b and 
b' in the qq component of the wave functions of the two onia. However, for the purpose of 
understanding the essentials of the behaviour of the elastic cross section, this complication 




(31) 
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Figure 9: The elastic and total cross sections for onium-onium scattering, as a function of 
rapidity, showing both the one-pomeron approximation and the fully unitarised results. 

will be ignored: \F\ will be used instead of |^|. The leading power behaviour of F is 1/r^, 
which is why for the total cross section, large impact parameters can contribute. For the 
elastic cross section, though, one is integrating the square of the amplitude, which will fall 
off as 1/r'^. So the integral is dominated by small r, where the unitarisation is strongest. 
Another aspect of the problem, is that the elastic cross section can't exceed half the total 
cross section (the unitarity condition), but the asymptotic power growth with Y of the elastic 
cross section (in the one pomeron approximation) is twice that of the total cross section. 
Therefore unitarisation must reduce the power growth of the elastic cross section by more 
than a factor of two. 

Figure shows the one pomeron and unitarised calculations for the elastic cross section, 
together with the analogous results for the total cross sections, for reference. As suggested, 
the unitarisation corrections to the elastic cross section are much larger than those to the total 
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Figure 10: The power dependence in rapidity of the total and elastic, 1-pomeron and uni- 
tarised cross sections, as a function of rapidity. 



cross section, and one clearly sees the rate of increase of the elastic cross section tailing off 
to take the same value as that of the total cross section. One also sees that at large rapidity, 
the fraction of the total cross section which comes from elastic scattering is relatively small 
(about 0.06), certainly well below the unitarity bound of 1/2. The reason for this is, again, 
that the total cross section is coming from a large region in impact parameter where the 
amplitude is small, and so does not contribute much to the elastic cross section. It is also 
worth briefly noting, in view of the discussion on eikonalisation in the previous section, that 
the rapidity where the one-pomeron total and elastic cross sections cross over (Y ~ 20), is 
the same as the point where the one and two pomeron contributions to the total cross section 
would cross if the eikonal approximation were valid. 



In figure the effective power dependence of the total and elastic cross sections is plotted 
against Y. Looking first at the 1 pomeron cross sections, one sees that they vary with Y, 
tending only relatively slowly to their asymptotic values of 0.47 and 0.94 respectively. This 
is a result of the logarithmic corrections which are 1/VY and 1/Y^ for the total and elastic 
cross sections. The power for the unitarised total cross section deviates only slowly from the 
1-pomeron power, and its variation with Y is also slow. The unitarised elastic cross section, on 
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the other hand, grows with a power which differs quite noticeably from the 1-pomeron elastic 
cross section, and the variation of the power with Y is also very different, decreasing rather 
than increasing from Y = 8 onwards, tending towards the same power as the unitarised total 
cross section at large rapidity. This significant decrease in the elastic power, if not masked 
by other effects, might be a clear signal for the onset of unitarity corrections. 
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Figure 11: The unitarised differential elastic cross section for different values of Y. 

Finally, for completeness, it is worth looking at the differential elastic cross section, fig- 
ure as a function of the square of the exchanged momentum, t. It is very strongly peaked 
at t = 0, because in impact parameter the amplitude dies off so slowly at large r. As Y in- 
creases, the region of interaction gets larger, which causes the elastic peak to become narrower 
in t. 

5 Conclusions 

By performing a Monte Carlo simulation of the dipole formulation of onium-onium elastic 
scattering, it has been possible to calculate single and multi-pomeron scattering amplitudes 
as a function of rapidity and impact parameter. The energy where one and two pomeron 
total amplitudes are comparable, has been determined to be y ~ log s/M^ ~ 14 (with 
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as = 0.18). However, the multiple pomeron series is not a very good way of calculating 
unitarity corrections. There are two reasons for this. Firstly, the multiple pomeron series 
diverges. This is a consequence of the dipole density in an onium wave function having 
exponential fluctuations. For multiple pomeron exchange, the enhancement due to multiple 
coupling to these very dense configurations more than compensates for their low probability. 
This necessitates the resummation of all numbers of pomeron exchange before averaging over 
onium configurations. Nevertheless, one might have expected that once, say, two pomeron 
exchange became larger than one pomeron exchange, the unitarised amplitude would quickly 
stop growing. It doesn't: unitarity corrections set in very slowly, with the correction to the 
effective power of the growth being only of the order of 25%, even at a rapidity of y = 19. 

The main reason for this is that the growth of the total multiple pomeron amplitude is 
due mostly to an increase in the interaction at central impact parameters. The unitarised 
amplitude, on the other hand effectively carries on growing through an increase in the area 
of interaction, even when it has already reached the unitarity bound at a given point in 
impact parameter. This though is still a simplification. One finds that the impact parameters 
contributing most to the total amplitude have low average amplitudes, yet significant unitarity 
corrections. In fact, the average amplitude there is dominated by rare configurations with 
large amplitudes and large unitarity corrections: the growth of the total amplitude comes 
through an increase in the probability of these rare configurations. 

Since the unitarisation is largest at small impact parameters, it is useful to find a process 
where small impact parameters matter more. Elastic scattering has this property because it 
depends on the square of the average amplitude so regions with the largest amplitudes (and 
strongest unitarisation) contribute most. One finds that the unitarisation effects start to 
become important as early as y = 8: the one-pomeron approximation to dlogUei/dt carries 
on increasing with Y, while the unitarised power starts to decrease. 

There are a number of difficulties when trying to relate these results to experimentally 
feasible measurements. Ideally one would want to examine a process such as that suggested 



by Mueller and Navelet |13|, where the BFKL evolution occurs between two objects with 
similar, large, transverse scales. However, the evidence currently for BFKL type behaviour 
is in results from HERA, both in the measurement of the rise of the structure function F2 at 
small X ||l^, |l5| (analogous to a total cross section) and more recently from the rise (roughly 
twice as fast) of the elastic electro-production of and J/t/; and elastic photo-production 
of J/ip |]16|, |l^, 19, ^] (analogous to elastic scattering). Complications arise in trying to 
understand the unitarity corrections in these types of processes, partly because the evolution is 
taking place over a range of scales, but also because of the non-perturbative effects associated 
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with the proton scale. In particular, these would prevent the large size fluctuations, which 
are so important for the continued growth of the BFKL total cross section beyond the point 
where unitarisation becomes significant. 
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Appendix A: The mean spatial distribution of dipoles. 



The discrepancy between the Monte Carlo results for the dipole distribution and the analytical 
prediction from Mueller is related to the details of one of the integrals involved in its 
evaluation. In determining eq. ([l^ ) one goes through the following stage: 



n{c,b,r,y)= [ ^H^{c,h,r)- exp 
7 zvr c 



2asNc 



IT 



X(l + 2iu)y 



(32) 



where 



2iu „2iu 



F^(c,b,r) = -^b-'^'c 



d2R 



|R - b/2||R + b/2|]i-2»'^[|R + r - c/2||R + r + c/2\]^+^^'' 

(33) 



The BFKL characteristic function, xi^)? is 



x(7)=V'(l)-^V'(l-7/2)-^V'(7/2), 



(34) 



where tp is the digamma function. First, cut out the regions R < Pb and |R + r| < Pc, where 
P is an arbitrary constant satisfying P ^ 1 and 4iyP <^ 1. The integral without these regions 
gives: 



2iiJ 



(—] 



-2iu 



(—] 



2iu 



(35) 



One can see, by taking the limit v ^ 0, that there should be no terms of order i/ within 
the square bracket, because the integral is similar to the one evaluated in for the virtual 
corrections to the dipole kernel. The result obtained in is equivalent to taking P = 1 and 
neglecting other contributions. Now consider one of the regions that was cut out. It can be 
integrated as follows: 



4^j^2jj-2iu^2tu .R<Ph RdRde 



r^+'ii'^ J {RVb-^ - 2bR cose + R^)^-^'^ 

4u^{PHc)^''' fR<P^ dRde 



(36) 

(•37) 

vr2 r2+4i^^ J Vb'' -2bRcose + R'' ^ ' 

The approximation giving the second line is valid because of the condition vP ^ 1 , so that in 
the range of R important for the integration, the power v has little effect on the value of the 
integrand. The angular integral can be performed to give a complete elliptic integral of the 



first kind, [4K(i?/6)/6] ioi R<b and [4K(6/i?)/i?] for > 6 H. Performing the R integral 
then gives 
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Adding this, and a similar term for the dipole of size c, to eq. p5|) gives 
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Following through the calculation, gives for the dipole density: 

b r y) - log ( exp[(ap - l)y - log^{16r^ /bc)/ky] 

This is much more central and agrees well with the Monte Carlo distributions. 
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Appendix B: Fluctuations in the localised dipole number 

The approach used, will be to study the equations for the moments of the dipole distribution 
within a region of radius p. The first step is to note that branching sequences that include 
dipoles of sizes larger than p will not contribute to high density fluctuations: the number 
dipoles of size p coming from a dipole of size R ^ p will be at most enhanced by a factor 
proportional to R/p (just the usual factor in front of the BFKL saddle point solution). How- 
ever dipoles will be spread out over a region of size R^, so that the density of dipoles will 
be suppressed by a factor p/R. Therefore in calculations one can leave out all dipoles larger 
than p. One can also neglect the transverse positions of the dipoles, since, if a dipole starts 
inside a region of size p, and the branching involves no dipoles larger than p, most of the child 
dipoles will still be inside the region. 

Define n^p^ (c, 6, y) to be the q^'^ moment of the number of dipoles of size c in a region of 
radius p originating from a parent of size b after an evolution through rapidity y. (The q^^ 
moment being (n(n — 1) ... (n — (7 + 1))). It will then approximately satisfy the following 
equation: 



d4''\e, b,y) _ asNc fP b^fh2 

dy ^ + 27r2 i bl,bl. 

The inhomogeneous term, Ip'^\ is 



4'\c,bo2,y) + 4''\c,bu,y)-4'Hc,boi,y) ■ (41) 



= ^ r'^^:c?nPicM2,y)n^r'ic,h2,y) 
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This is obtained from the generating functional equation in Q, with an upper cutoff placed 
on the dipole size in the b2 integration. 

The first approximation to make relies on the large Y limit. At large Y, with c and b not 
too different, one has for the total average number of dipoles: 



n{c,b,Y) oc -e("^-i)'^ 



(43) 



The important point is that the b dependence (which shows a power behaviour) factors out 
from the c and Y dependences. This will apply to the moments of the localised dipole 
distribution as well. It will further be useful to make the assumption (reasonably justifiable) 
that the b dependences of the moments are also power behaviours. One can then parameterise 
each moment by two constants and Vq, which should depend only weakly on b: 



n^;^'{c,b,Y)=eqq\A\c/p) 



(44) 



It is not necessary to know the exact c dependence, which is why it is left as A{c). With this 
parameterisation, eq. ( pl| ) can be written as: 



^P'^ i=l 



(45) 



The following integrals have been defined: 
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For 1 < < 2 (as should be the case for all q), ky and lyyt will depend only weakly on b/r 
(except for In which has a logarithmic dependence on b/r — this is discussed later). 

The main difference between the two terms is that the l^iy' integral is dominated by large 
62 ~ p, while the ku integral, because of the condition given above on z^, the integral is 
dominated by small 62 ~ b. This is the origin of the different power behaviours from the two 
integrals. 

To solve eq. ( ^5| ) one uses the fact that < 2, therefore for large p/b the inhomogeneous 
term drops out, leaving: 



j{av - 1) 



asNc 



ki/„ 



(48) 
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Through the v dependence in ky, this fixes Vq. For moderately large q, ky can be approximated 
as 

^ (49) 

giving 

1 



- - 2, log 2 

Note that this is just the solution to the equation x(^) = iii the limit ^ 2, with 

X defined as in eq. ([3^). The other region which can be studied to help solve eq. (^) is 
that of h close to p. For the case of v close to 2, i.e. for large q, at large 6, the integral 
for ku is much smaller than its asymptotic value (it only reaches its asymptotic value for 
p/h ~ exp[l/(2 — y)]), while that for l^j^i should be much closer to its asymptotic value (the 
integrand is approximately proportional to 62)- Therefore the ky term can be neglected, giving 



1 1 



^ 41og2g 
A reasonable approximation for lu^u2 is: 



'^^eieq-ily^l,^_^ (51) 



By iterating eq. (51), one can determine all the moments of the fluctuations. Its form is such 



that for large q, Cg should be behave as follows: 

eq ~ BC"? (53) 

Because of the complicated behaviour of the coefficients of the iteration for low values of q, 
it is not possible to provide an analytical approximation for B and C, though they can be 
determined numerically. This then gives an exponential distribution for the fluctuations in 
dipole number in a restricted region, e~"/ ^ where 

/j{c/r) = A(c/r)Ce("^-i)^ (54) 

It is difficult to compare this directly with the results from the Monte Carlo simulation, 
because A(c/r) contains various logarithmic factors which are difficult to estimate accurately. 
However the mean number of dipoles of a given size in a certain region is something that can 
be obtained from the Monte Carlo distribution and so A{c/p) can be factored out: 
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In evaluating C from the recursion relations there are some complications because the first 
relation contains In which varies logarithmically with h/ p (this is the reason why the expres- 
sion given in eq. ( |52|) diverges). Given that the values of Cq are determined by the region of 
h close to p, it is then natural to use a value of In corresponding also that region of b. For 
h = 0.5p, III = 1-9 (determined numerically), while for b = 0.25/3, In = 2.7. This range will 
be used when calculating the eg to gauge the uncertainty in the calculation. 

Another uncertainty arises because at the rapidities available to the Monte Carlo simu- 
lation, the logarithmic prefactors affect the effective power growth of the distribution. To 
gauge this effect, one can replace {ap — 1) with a measured effective power and look at the 
difference that this makes to C. The final range of values obtained for C is from 0.34 to 0.75. 
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Figure 12: The parameter C of eq. ([55|), as determined by Monte Carlo simulation for two 
values of y. 




Figure ^ shows the value of Cmc = A*(c/p)/[2np(c, p/2, y)] determined from the Monte 
Carlo simulation for two values of Y and a range of c/ p. The values for C from the iteration 
of eq. (51) are consistent with the Monte Carlo results. However the Monte Carlo results do 
depend slightly on Y and on c/p. For the variation with Y at c ^ p, Cmc decreases as Y 
increases. This could be because at larger Y, the effective power is larger. A larger effective 
power leads to a reduced value for eg/cq-i (this can be seen by making some simple analytic 
approximations, and is also obtained in the numerical iterations). A second effect is that 
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at smaller c/p, the presence, non-asymptotically, of factors of the form exp[— log^(c//9)/A;y] 
leads to an increase in the effective power with decreasing c, giving a decreasing C. At larger 
Y this effect would be smaller, so, at larger Y, C should decrease more slowly with decreasing 
c. Both of these observations are qualitatively in accord with the Monte Carlo results. 
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